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Abstract. We present the results of our analysis of a high resolution {R ~ 30000) infrared spectrum of V4334 Sgr 
(Sakurai's Object) around the first overtone CO bands, obtained in 1998 July. The ^^CO and ^^CO bands are 
well-resolved, and we compute synthetic hydrogen-deficient model atmosphere spectra to determine the ^^C/^^C 
ratio. We find ^^C/^''C ~ 4 ± 1, consistent with the interpretation of V4334 Sgr as an object that has undergone 
a very late thermal pulse. 



Key words, stars: individual: V4334 Sqr 
stars: abundances 



stars: individual: Sakurai's Object - stars: atmospheric parameters 



1. Introduction 

V4334 Sgr (Sakurai's Object) is a low mass (~ 0.6- 
O.8M0; Jlcrwie 2001) star undergoing a very late thermal 
pulse (VLTP). After traversing the 'knee' on the post- 
asymptotic giant branch (AGB) evolutionary track, the 
star experiences a final He flash as it evolves down the 
white d warf coohng track: it becomes a 'born-again ' AGB 
star (see lHerwidf200lULawlor fc MacDonaldll2003(l . 

A key factor in confirming this scenario is the determi- 
nation of elemen tal abundances and iso topic ratios in the 
surface material. lAsplund et alJ l)l999(l have used optical 
spectroscopy to determine the abundance of a range of 
species; they found V4334 Sgr to be hydrogen-deficient, 
with other elemental abundances close to those of the 
R CrB stars; R CrB stars. They also used the C2 (1-0) 
and (0-1) Swan bands at 4740 and 5635 A respectively to 
determine that the ^'^C/^^C ratio was somewhere in the 
range 1.5-5, although they did not give an actual value. 

IPavlenko fc Geballd l|2002f) modelled the spectral en- 
ergy distribution of V4334 Sgr in the near -IR (AA 1 - 2.5 
um). They broadly confirmed the results of lAsplund et alJ 
l|l999l) . but foun d some evidences of h igher C abundance, 
by a factor ~ 2.'Pavlenk o fc Geballel also noted the pres- 
ence of hot dust longward of ~ 2 /im. 

The ^^C/^^C ratio is of fundamental importance in 
understanding the VLTP process, as it is a refiection of 
the relative importance of ^^C, ^^C processing and dredge- 
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up. In this Letter we use echelle observations of the first 
overtone vibration-rotation features of CO, obtained in 
1998 July while the star was still visible, to determine the 
i^C/i^C ratio in V4334 Sgr. 



2. Observation 

The data were obtained with th e facility spectrome- 
ter CGS4 l|Moimtain et all Il99nl) at the 3.8m United 
Kingdom Infrared Telescope on UT 1998 July 7. CGS4 
was configured with a 256 x 256 InSb array, its 31 Imm^^ 
echelle and a 0.45" slit, which provided a resolving power 
of ~ 30, 000. The spectrum was sampled every 0.4 resolu- 
tion element. Five echelle settings were used to continu- 
ously cover the wavelength range 2.32-2.38 /^m, an inter- 
val containing numerous ^^CO and ^■^CO lines, including 
three bandheads of ^^CO ((5u = 3 ^ 1, 4 ^ 2, 5 ^ 3) 
and two of "CO {Sv 2 ^ 0,3 1). Each of the five 
spectral segments was observed for a total of 8 minutes, 
with individual data points being observed for 4 minutes. 
A spectrum of the A2V star HR6378 was obtained just 
prior to each spectrum of V4334 Sgr and at similar air- 
mass. The raw spectra of V4334 Sgr and IIR6378 were ex- 
tracted, reduced, and ratioed using standard techniques. 
Wavelength calibration was achieved using some of the 
multitude of telluric absorption lines in the spectra of the 
calibration star, and is accurate to ±3kms^^; the wave- 
length scale is in vacuum. The ratioed spectral segments 
slightly overlap one another and were scaled and adjoined 
to provide the final spectrum used in the analysis. 
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3. Procedure 

3.1. Spectral synthesis 

We have used the technique of synthetic spectra to carry 
out our analysis of the IR spectrum of V4334 Sgr. We 
computed the spectra in the classical framework, assum- 
ing LTE, plane-parallel media, and no sinks and sources of 
energy in the atmosphere. The transfer of energy is pro- 
vided by the radiation field and by convection. Strictly 
speaking, none of these assumptions is completely valid 
in the atmosphere of V4334 Sgr because, in all probabil- 
ity, we see only the pseudophotosphere, i.e. the outermost 
part of the expanding envelope. However, we suggest that 
none of these assumptions is of crucial importance for the 
spectrum formation processes. 

A grid of LTE synthetic spectra was computed 
for 12C/13C ratios of 1, 2, 3, 4, 5, 6, 7, 8, 9, 
10, model atmospher e param eters Tefi,log5 = 5250,0.0 
I Pavlenko fc DiirbeckL l2002t) . in the wavelength range 
2.31-2.39 /im and with wavelength step 0.00002 ^m, us- 
ing the WITA612 co de llPavlenkolll997l). W e used abun- 
dances as given by lAsplund et al.l l)l999|) . scaled such 
that = ^^'^ ■^itli O abundances of -2.52, - 

2.02 respectively, and two abundances of carbon, namely 
log7V(C) = -1.05 and log N{C) = -1.62; the first 
case corresponds to a carbon- rich model atmosphere (see 
lAsplund et al.lll99^ for a discussion of the carbon prob- 
lem). As the determination of t he ^^C/^'^C ratio in st ellar 
atmospheres depends on Vt fsee IPavlenko et a'l1l2003l) . we 
determined ^^C/^^C for two values of the microturbulcnt 
velocity, Vt= 3 and 6kms~^. 

A se t of co ntinuum opacity sources from ATLAS9 
llKuruc^ Il993^ was used. Additionally, all the major 
sources of opacity in our spectral region of interest were 
taken into account: 

(i) bo und-free absorption by C I, Nij O i atoms computed 
by IPavlenko fc Zhukovska' ( 2003) for the TOPBASE 
cross-sections (,Seaton et al. ,1992): 

(ii) 12C160 , lines from lCxOorvitd] lll994l) : 

(iii) iZQi ^C, i^C^C, " C13C, 12C"N and "c^n lines 
from lKuruc3 l)l993|) : 

(iv) atomic lines from VALD l|Kupka et alill999|) . 

We adopt the Voigt profile H{a,v) for the shape of 
each line; damping constants a = (72 -I- 74 -I- 76)/(4 x tt x 
Avn) were computed using data f rom various databases 
ijKurucj. I1993I: iKupka et allll999tl. o r computed follow- 
ing an approximation of' IUns51dri)l955(l . The contributions 
of the various molecular and atomic species to the total 
opacity in the 2.31-2.39 /im spectral region are shown in 
Fig. n from which we see that absorption by CO dom- 
inates. Finally, beari ng in rnind th e high luminosity of 
V4334 Sgr (- 10"^ L^: lHerWi3l200lh . we adopt logg = 0. 

3.2. The effect of circumstellar dust 

IPavlenko fc aeballel 1)2002') have shown that the 2 /im 
spectral region of V4334 Sgr is affected by emission from 
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Fig. 1. Theoretical spectra of different species in the spec- 
tra region convolved by gaussian of half- width 5x10^^ /im. 
Synthetic spectra were computed for ^^0/^^0=1 

the hot dust known to envelop V4334 Sgr at this time, 
and which most likely fo rmed at some 20 stellar radii from 
the stellar atmosphere (iKippeil Il999t iTvne et all l2002^ . 
The dust shell first appeared in 1997 March and, by 1998 
August, its presence was wel l estab lished photometrically 
tee rber et al ] '1999': 'K ipped Il999l) and spectroscopically 
()Tvne et al.l I200Q. .200^ 

The dust emission veils the CO absorption lines, and to 
estimate its effect we suppose that the dust emission adds 
some continuum fiux in our spectral region, effectively re- 
ducing the residual fluxes r^, — F}^'^ /F^ (and hence equiv- 
alent widths) of the absorption lines; here Fl'^^ and F^ are 
fluxes in the lines and continuum respectively. We describe 
this by the factor 

= (rO+eo)/(l + eo), 

where and eo are respectively the residual flux com- 
puted for the 'dust-free' case, and the relative contribu- 
tion of the dust emission; the parameter eq is determined 
from comparison of computed and observed spectra. 

3.3. Fitting procedure 

To determine the best flt parameters, we compare the ob- 
served residua l fluxes r„ w i th com put ed values following 
the scheme of Ijones et al.l l|2002|) and IPavlenko fc Jonea 
( 2002) . We let 

= j rly^ G{x-y)*dy, 

where rf, and G{x — y) are respectively the residual fluxes 
computed by WITA612 and the broadening profile; we 
adopt a gaussian for the latter. We then find the minima 
of the 3D function 

S{h, /h, /g) = E (1 - X r^y"Vr°'^^)' , 
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Fig. 2. a) Dependence of computed spectra on ^^C/^'^C for values of Toff, \ogg and eo indicated, b) Dependence of 
the parameter S on ^^C/^^C for values of logiV(C) and Vt indicated. 



Table 1. Atmospheric and dust parameters for V4334 Sgr 
in 1998 July. 



logiV(C) 
Vt (kms-i) 


3 


1.05 
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1.62 
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Shift (10~*/.im) between 










observed and computed spectra 


8.6 


8.6 


8.6 


8.6 


Normalisaton factor /h 


0.96 


0.96 


0.96 


0.96 


FWHM (10-* pm) 


1.56 


1.4 


1.8 


1.4 


eo 


1.4 


2.0 


1.8 


2.4 


12^/13^ 


4 


3 


(5,8) 


5 



where fs, /g, /g are respectively the wavelength shift, nor- 
malisation factor, and half- width of G{x — y). The param- 
eters /s, /h and /g are determined by the minimization 
procedure for every computed spectrum. Then, from the 
grid of the better solutions for a given ^^C/^'^C and eg, we 
choose the best-fitting solution. 

4. Results 

On July 7 the Earth's orbital motion redshifted the 
spectrum of V4334 Sgr by 8kms"^ Thus the shift of 
-fll2kms^^ that we needed to apply means that the 
heliocentric velocity of the photosperic CO was -1-104 ± 
3.2 km s"^, where the uncertainty includes ±3kms~^ 
for the wavelength scale, and ±lkms^^ in our shift. 
While this is close to the known h eliocentric velocity of 
V4334 Sgr, -1-115 km s'^ l)Diirbeck fc Benettl 1996), it 
may also imply that the photospheric layers in which the 
CO was located were expanding at ~ lOkms"^. 

Our results show that ^^C/^^C depends on several pa- 
rameters. Fig. |3l shows the dependence of the synthetic 
spectra on the ^^C/^^C ratio for eo = 2.0, while Fig. |2b 
shows the dependence of the best-fit parameter S on the 
12q^13q Yaiio for a variety of logiV(C) and Vt values. Both 
the ^^C^^O and ^'^C^^O bands are strong in the spectra 
due to the high carbon abundance; the computed central 



intensities of the strongest lines are ^ 0.4 of the normal- 
ized continuum. 

Using our fitting procedure we found a clear depen- 
dence of the spectra on ^^C/^'^C, Vt and eo. The best- 
fit synthetic spectra for our two assumed C abundances 
are shown in Figs.|31i,b; the corresponding parameters are 
given in Tabled Note that there are two possible solutions 
for [logiV(C), Vt] = [-1.62, 3] (see Fig.l^s). Fig.|ii,b shows 
the 2D-dependence of the fit parameter S on ^^C/^^C and 
eo- We effectively found two families of solutions, differing 
by their ^^C/^^C ratios and eo values. An alternative in- 
terpretation is that there is an uncertainty in our results 
dependent on the uncertainties in the input parameters. 

In general, a mictoturbulent velocity Vt= 6kms-^ is 
most likely for a supergiant atmosphere and our results for 
Vt= 3kms-^, presented in Tabled are included primarily 
to illustrate the Vt-dependence of our results. We see from 
FigsEt andU and Tabled that a i^C/isC ratio of 4 ± 1 
is consistent with the data. 

5. Discussion 

One of the major challenges posed by V4334 Sgr has been 
the extrem e rapidi t y with which it has evolved following 
the VLTP. E erwid 1)200 1m has addressed this problem by 
supposing that the element mixing efficiency in the He 
flash convection zone during the VLTP is much smaller 
than predicted by mixing length theory. Herwig defines the 
parameter / = £'mlt/-Dcm, where I?mlt is the diffusion 
coefficient for convecti ve element mi xing and T'cm that 
for composition mixing . FHerwid (fioOl*) finds that the rapid 
evolution of V4334 Sgr is accounted for by / ~ 10 — 100, 
and that, for this value of /, the surface -'^^C/^^C ratio for 
the / = 30 VLTP model is 5. This is clearly very much 
in line with the results p resent ed here. 

IPavlenko & Ceballel l(2002l) found that dust made a 
major contribution to the emission longward of 2 /im in 
1998, and indeed was present as early as 1997. The value 
of the parameter eo gives an indication of the level of 
the dust emission at 2.3 /im. From Table ^ we see that 
eo — 2 ± 0.2, indicating that the dust flux at 2.3 /um is 
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Fig. 3. a) Best fit to 1998 July spectrum found by tlie minimisation procedure outlined in the text, for logiV(C) 
-1.05. b) As (a), but for logA/'CC) = -1.62 
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Fig. 4. a) Dependence of the fit parameter S on ^"^C/^^C and eq for a grid of the best fits found by minimisation 
procedure for logA^(C) = -1.05. b) As (a), but for logiV(C) = -1.62. 



~ 9.4 X 10^^^ W m^^ ^m^^. This is much higher than in 
later spectra, when the star had become completely ob- 
scured and the dust ein ission dominated the overall spec- 
trum l|Tvne et alil2002l) . At later times ( > 1999) the dust 
envelope expanded and cooled and the peak of its spec- 
trum shifted to longer wavelengths. 

6. Conclusion 

We have carried out a spectral synthesis to determine the 
i2(-;^i3(j j.g^^jQ V4334 Sgr. Given a range of C abun- 
dance, we determine that ^^C/^^C ~ 4±1, consistent with 
the VLTP interpretation of V4334 Sgr. The relative dust 
contribution at 2.3 /im is eo — 2 ± 0.2. Our value of the 
12(-i^13q patio is quite well constrained, and will doubtless 
be useful in fine-tuning the VLTP scenario. 
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